The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding.

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

An improved reference of the grapevine genome reasserts the origin of the PN40024 highly homozygous genotype.

The genome sequence of the diploid and highly homozygous Vitis vinifera genotype PN40024 serves as the reference for many grapevine studies. Despite several improvements to the PN40024 genome assembly, its current version PN12X.v2 is quite fragmented and only represents the haploid state of the genome with mixed haplotypes. In fact, being nearly homozygous, this genome contains several heterozygous regions that are yet to be resolved. Taking the opportunity of improvements that long-read sequencing technologies offer to fully discriminate haplotype sequences, an improved version of the reference, called PN40024.v4, was generated. Through incorporating long genomic sequencing reads to the assembly, the continuity of the 12X.v2 scaffolds was highly increased with a total number decreasing from 2,059 to 640 and a reduction in N bases of 88%. Additionally, the full alternative haplotype sequence was built for the first time, the chromosome anchoring was improved and the number of unplaced scaffolds was reduced by half. To obtain a high-quality gene annotation that outperforms previous versions, a liftover approach was complemented with an optimized annotation workflow for Vitis. Integration of the gene reference catalogue and its manual curation have also assisted in improving the annotation, while defining the most reliable estimation of 35,230 genes to date. Finally, we demonstrated that PN40024 resulted from 9 selfings of cv. "Helfensteiner" (cross of cv. "Pinot noir" and "Schiava grossa") instead of a single "Pinot noir". These advances will help maintain the PN40024 genome as a gold-standard reference, also contributing toward the eventual elaboration of the grapevine pangenome.

Characterization of Grapevine Fanleaf Virus Isolates in 'Chardonnay' Vines Exhibiting Severe and Mild Symptoms in Two Vineyards.

Fanleaf degeneration is a complex viral disease of Vitis spp. that detrimentally impacts fruit yield and reduces the productive lifespan of most vineyards worldwide. In France, its main causal agent is grapevine fanleaf virus (GFLV). In the past, field experiments were conducted to explore cross-protection as a management strategy of fanleaf degeneration, but results were unsatisfactory because the mild virus strain negatively impacted fruit yield. In order to select new mild GFLV isolates, we examined two old 'Chardonnay' parcels harbouring vines with distinct phenotypes. Symptoms and agronomic performances were monitored over the four-year study on 21 individual vines that were classified into three categories: asymptomatic GFLV-free vines, GFLV-infected vines severely diseased and GFLV-infected vines displaying mild symptoms. The complete coding genomic sequences of GFLV isolates in infected vines was determined by high-throughput sequencing. Most grapevines were infected with multiple genetically divergent variants. While no specific molecular features were apparent for GFLV isolates from vines displaying mild symptoms, a genetic differentiation of GFLV populations depending on the vineyard parcel was observed. The mild symptomatic grapevines identified during this study were established in a greenhouse to recover GFLV variants of potential interest for cross-protection studies.

Comparative Plant Transcriptome Profiling of Arabidopsis thaliana Col-0 and Camelina sativa var. Celine Infested with Myzus persicae Aphids Acquiring Circulative and Noncirculative Viruses Reveals Virus- and Plant-Specific Alterations Relevant to Aphid Feeding Behavior and Transmission.

Evidence is accumulating that plant viruses alter host plant traits in ways that modify their insect vectors' behavior. These alterations often enhance virus transmission, which has led to the hypothesis that these effects are manipulations caused by viral adaptation. However, we lack a mechanistic understanding of the genetic basis of these indirect, plant-mediated effects on vectors, their dependence on the plant host, and their relation to the mode of virus transmission. Transcriptome profiling of Arabidopsis thaliana and Camelina sativa plants infected with turnip yellows virus (TuYV) or cauliflower mosaic virus (CaMV) and infested with the common aphid vector Myzus persicae revealed strong virus- and host-specific differences in gene expression patterns. CaMV infection caused more severe effects on the phenotype of both plant hosts than did TuYV infection, and the severity of symptoms correlated strongly with the proportion of differentially expressed genes, especially photosynthesis genes. Accordingly, CaMV infection modified aphid behavior and fecundity more strongly than did infection with TuYV. Overall, infection with CaMV, relying on the noncirculative transmission mode, tends to have effects on metabolic pathways, with strong potential implications for insect vector-plant host interactions (e.g., photosynthesis, jasmonic acid, ethylene, and glucosinolate biosynthetic processes), while TuYV, using the circulative transmission mode, alters these pathways only weakly. These virus-induced deregulations of genes that are related to plant physiology and defense responses might impact both aphid probing and feeding behavior on infected host plants, with potentially distinct effects on virus transmission. IMPORTANCE Plant viruses change the phenotype of their plant hosts. Some of the changes impact interactions of the plant with insects that feed on the plants and transmit these viruses. These modifications may result in better virus transmission. We examine here the transcriptomes of two plant species infected with two viruses with different transmission modes to work out whether there are plant species-specific and transmission mode-specific transcriptome changes. Our results show that both are the case.

Comparative Metabolomic Analysis of Four Fabaceae and Relationship to In Vitro Nematicidal Activity against Xiphinema index.

The grapevine fanleaf virus (GFLV), responsible for fanleaf degeneration, is spread in vineyards by the soil nematode Xiphinema index. Nematicide molecules were used to limit the spread of the disease until they were banned due to negative environmental impacts. Therefore, there is a growing interest in alternative methods, including plant-derived products with antagonistic effects to X. index. In this work, we evaluated the nematicidal potential of the aerial parts and roots of four Fabaceae: sainfoin (Onobrychis viciifolia), birdsfoot trefoil (Lotus corniculatus), sweet clover (Melilotus albus), and red clover (Trifolium pratense), as well as that of sainfoin-based commercial pellets. For all tested plants, either aerial or root parts, or both of them, exhibited a nematicidal effect on X. index in vitro, pellets being as effective as freshly harvested plants. Comparative metabolomic analyses did not reveal molecules or molecule families specifically associated with antagonistic properties toward X. index, suggesting that the nematicidal effect is the result of a combination of different molecules rather than associated with a single compound. Finally, scanning electron microscope observations did not reveal the visible impact of O. viciifolia extract on X. index cuticle, suggesting that alteration of the cuticle may not be the primary cause of their nematicidal effect.

Differential gene expression in aphids following virus acquisition from plants or from an artificial medium.

BACKGROUND: Poleroviruses, such as turnip yellows virus (TuYV), are plant viruses strictly transmitted by aphids in a persistent and circulative manner. Acquisition of either virus particles or plant material altered by virus infection is expected to induce gene expression deregulation in aphids which may ultimately alter their behavior. RESULTS: By conducting an RNA-Seq analysis on viruliferous aphids fed either on TuYV-infected plants or on an artificial medium containing purified virus particles, we identified several hundreds of genes deregulated in Myzus persicae, despite non-replication of the virus in the vector. Only a few genes linked to receptor activities and/or vesicular transport were common between the two modes of acquisition with, however, a low level of deregulation. Behavioral studies on aphids after virus acquisition showed that M. persicae locomotion behavior was affected by feeding on TuYV-infected plants, but not by feeding on the artificial medium containing the purified virus particles. Consistent with this, genes potentially involved in aphid behavior were deregulated in aphids fed on infected plants, but not on the artificial medium. CONCLUSIONS: These data show that TuYV particles acquisition alone is associated with a moderate deregulation of a few genes, while higher gene deregulation is associated with aphid ingestion of phloem from TuYV-infected plants. Our data are also in favor of a major role of infected plant components on aphid behavior.

Severe Stunting Symptoms upon Nepovirus Infection Are Reminiscent of a Chronic Hypersensitive-like Response in a Perennial Woody Fruit Crop.

Virus infection of plants can result in various degrees of detrimental impacts and disparate symptom types and severities. Although great strides have been made in our understanding of the virus-host interactions in herbaceous model plants, the mechanisms underlying symptom development are poorly understood in perennial fruit crops. Grapevine fanleaf virus (GFLV) causes variable symptoms in most vineyards worldwide. To better understand GFLV-grapevine interactions in relation to symptom development, field and greenhouse trials were conducted with a grapevine genotype that exhibits distinct symptoms in response to a severe and a mild strain of GFLV. After validation of the infection status of the experimental vines by high-throughput sequencing, the transcriptomic and metabolomic profiles in plants infected with the two viral strains were tested and compared by RNA-Seq and LC-MS, respectively, in the differentiating grapevine genotype. In vines infected with the severe GFLV strain, 1023 genes, among which some are implicated in the regulation of the hypersensitive-type response, were specifically deregulated, and a higher accumulation of resveratrol and phytohormones was observed. Interestingly, some experimental vines restricted the virus to the rootstock and remained symptomless. Our results suggest that GFLV induces a strain- and cultivar-specific defense reaction similar to a hypersensitive reaction. This type of defense leads to a severe stunting phenotype in some grapevines, whereas others are resistant. This work is the first evidence of a hypersensitive-like reaction in grapevine during virus infection.

Alternative splicing regulation appears to play a crucial role in grape berry development and is also potentially involved in adaptation responses to the environment.

BACKGROUND: Alternative splicing (AS) produces transcript variants playing potential roles in proteome diversification and gene expression regulation. AS modulation is thus essential to respond to developmental and environmental stimuli. In grapevine, a better understanding of berry development is crucial for implementing breeding and viticultural strategies allowing adaptation to climate changes. Although profound changes in gene transcription have been shown to occur in the course of berry ripening, no detailed study on splicing modifications during this period has been published so far. We report here on the regulation of gene AS in developing berries of two grapevine (Vitis vinifera L.) varieties, Gewurztraminer (Gw) and Riesling (Ri), showing distinctive phenotypic characteristics. Using the software rMATS, the transcriptomes of berries at four developmental steps, from the green stage to mid-ripening, were analysed in pairwise comparisons between stages and varieties. RESULTS: A total of 305 differential AS (DAS) events, affecting 258 genes, were identified. Interestingly, 22% of these AS events had not been reported before. Among the 80 genes that underwent the most significant variations during ripening, 22 showed a similar splicing profile in Gw and Ri, which suggests their involvement in berry development. Conversely, 23 genes were subjected to splicing regulation in only one variety. In addition, the ratios of alternative isoforms were different in Gw and Ri for 35 other genes, without any change during ripening. This last result indicates substantial AS differences between the two varieties. Remarkably, 8 AS events were specific to one variety, due to the lack of a splice site in the other variety. Furthermore, the transcription rates of the genes affected by stage-dependent splicing regulation were mostly unchanged, identifying AS modulation as an independent way of shaping the transcriptome. CONCLUSIONS: The analysis of AS profiles in grapevine varieties with contrasting phenotypes revealed some similarity in the regulation of several genes with developmental functions, suggesting their involvement in berry ripening. Additionally, many splicing differences were discovered between the two varieties, that could be linked to phenotypic specificities and distinct adaptive capacities. Together, these findings open perspectives for a better understanding of berry development and for the selection of grapevine genotypes adapted to climate change.

Evidence of differential spreading events of grapevine pinot Gris virus in Italy using datamining as a tool.

Since its identification in 2003, grapevine Pinot gris virus (GPGV, Trichovirus) has now been detected in most grape-growing countries. So far, little is known about the epidemiology of this newly emerging virus. In this work, we used datamining as a tool to monitor in-silico the sanitary status of three vineyards in Italy. All data used in the study were recovered from a work that was already published and for which data were publicly available as SRA (Sequence Read Archive, NCBI) files. While incomplete, knowledge gathered from this work was still important, with evidence of differential accumulation of the virus in grapevine according to year, location, and variety-rootstock association. Additional data regarding GPGV genetic diversity were collected. Some advantages and pitfalls of datamining are discussed.

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication.

BACKGROUND: A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized. RESULTS: We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined. CONCLUSIONS: We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

Grapevine virus T diversity as revealed by full-length genome sequences assembled from high-throughput sequence data.

RNASeq or double-stranded RNA based approaches allowed the reconstruction of a total of 9 full-length or near full-length genomes of the recently discovered grapevine virus T (GVT). In addition, datamining of publicly available grapevine RNASeq transcriptome data allowed the reconstruction of a further 14 GVT genomes from five grapevine sources. Together with four GVT sequences available in Genbank, these novel sequences were used to analyse GVT diversity. GVT shows a very limited amount of indels variation but a high level of nucleotide and aminoacid polymorphism. This level is comparable to that shown in the closely related grapevine rupestris stem pitting-associated virus (GRSPaV). Further analyses showed that GVT mostly evolves under conservative selection pressure and that recombination has contributed to its evolutionary history. Phylogenetic analyses allow to identify at least seven clearly separated groups of GVT isolates. Analysis of the only reported PCR GVT-specific detection primer pair indicates that it is likely to fail to amplify some GVT isolates. Taken together these results point at the distinctiveness of GVT but also at the many points it shares with GRSPaV. They constitute the first pan-genomic analysis of the diversity of this novel virus.

Assessment of a Targeted Gene Panel for Identification of Genes Associated With Movement Disorders.

Importance: Movement disorders are characterized by a marked genotypic and phenotypic heterogeneity, complicating diagnostic work in clinical practice and molecular diagnosis. Objective: To develop and evaluate a targeted sequencing approach using a customized panel of genes involved in movement disorders. Design, Setting and Participants: We selected 127 genes associated with movement disorders to create a customized enrichment in solution capture array. Targeted high-coverage sequencing was applied to DNA samples taken from 378 eligible patients at 1 Luxembourgian, 1 Algerian, and 25 French tertiary movement disorder centers between September 2014 and July 2016. Patients were suspected of having inherited movement disorders because of early onset, family history, and/or complex phenotypes. They were divided in 5 main movement disorder groups: parkinsonism, dystonia, chorea, paroxysmal movement disorder, and myoclonus. To compare approaches, 23 additional patients suspected of having inherited cerebellar ataxia were included, on whom whole-exome sequencing (WES) was done. Data analysis occurred from November 2015 to October 2016. Main Outcomes and Measures: Percentages of individuals with positive diagnosis, variants of unknown significance, and negative cases; mutational frequencies and clinical phenotyping of genes associated with movement disorders. Results: Of the 378 patients (of whom 208 were male [55.0%]), and with a median (range) age at disease onset of 31 (0-84) years, probable pathogenic variants were identified in 83 cases (22.0%): 46 patients with parkinsonism (55% of 83 patients), 21 patients (25.3%) with dystonia, 7 patients (8.4%) with chorea, 7 patients (8.4%) with paroxysmal movement disorders, and 2 patients (2.4%) with myoclonus as the predominant phenotype. Some genes were mutated in several cases in the cohort. Patients with pathogenic variants were significantly younger (median age, 27 years; interquartile range [IQR], 5-36 years]) than the patients without diagnosis (median age, 35 years; IQR, 15-46 years; P = .04). Diagnostic yield was significantly lower in patients with dystonia (21 of 135; 15.6%; P = .03) than in the overall cohort. Unexpected genotype-phenotype correlations in patients with pathogenic variants deviating from the classic phenotype were highlighted, and 49 novel probable pathogenic variants were identified. The WES analysis of the cohort of 23 patients with cerebellar ataxia led to an overall diagnostic yield of 26%, similar to panel analysis but at a cost 6 to 7 times greater. Conclusions and Relevance: High-coverage sequencing panel for the delineation of genes associated with movement disorders was efficient and provided a cost-effective diagnostic alternative to whole-exome and whole-genome sequencing.

Annotation, classification, genomic organization and expression of the Vitis vinifera CYPome.

Cytochromes P450 are enzymes that participate in a wide range of functions in plants, from hormonal signaling and biosynthesis of structural polymers, to defense or communication with other organisms. They represent one of the largest gene/protein families in the plant kingdom. The manual annotation of cytochrome P450 genes in the genome of Vitis vinifera PN40024 revealed 579 P450 sequences, including 279 complete genes. Most of the P450 sequences in grapevine genome are organized in physical clusters, resulting from tandem or segmental duplications. Although most of these clusters are small (2 to 35, median = 3), some P450 families, such as CYP76 and CYP82, underwent multiple duplications and form large clusters of homologous sequences. Analysis of gene expression revealed highly specific expression patterns, which are often the same within the genes in large physical clusters. Some of these genes are induced upon biotic stress, which points to their role in plant defense, whereas others are specifically activated during grape berry ripening and might be responsible for the production of berry-specific metabolites, such as aroma compounds. Our work provides an exhaustive and robust annotation including clear identification, structural organization, evolutionary dynamics and expression patterns for the grapevine cytochrome P450 families, paving the way to efficient functional characterization of genes involved in grapevine defense pathways and aroma biosynthesis.

Combinatorial DNA methylation codes at repetitive elements.

DNA methylation is an essential epigenetic modification, present in both unique DNA sequences and repetitive elements, but its exact function in repetitive elements remains obscure. Here, we describe the genome-wide comparative analysis of the 5mC, 5hmC, 5fC, and 5caC profiles of repetitive elements in mouse embryonic fibroblasts and mouse embryonic stem cells. We provide evidence for distinct and highly specific DNA methylation/oxidation patterns of the repetitive elements in both cell types, which mainly affect CA repeats and evolutionarily conserved mouse-specific transposable elements including IAP-LTRs, SINEs B1m/B2m, and L1Md-LINEs. DNA methylation controls the expression of these retroelements, which are clustered at specific locations in the mouse genome. We show that TDG is implicated in the regulation of their unique DNA methylation/oxidation signatures and their dynamics. Our data suggest the existence of a novel epigenetic code for the most recently acquired evolutionarily conserved repeats that could play a major role in cell differentiation.

TRF2 is recruited to the pre-initiation complex as a testis-specific subunit of TFIIA/ALF to promote haploid cell gene expression.

Mammalian genomes encode two genes related to the TATA-box binding protein (TBP), TBP-related factors 2 and 3 (TRF2 and TRF3). Male Trf2(-/-) mice are sterile and characterized by arrested spermatogenesis at the transition from late haploid spermatids to early elongating spermatids. Despite this characterization, the molecular function of murine Trf2 remains poorly characterized and no direct evidence exists to show that it acts as a bona fide chromatin-bound transcription factor. We show here that Trf2 forms a stable complex with TFIIA or the testis expressed paralogue ALF chaperoned in the cytoplasm by heat shock proteins. We demonstrate for the first time that Trf2 is recruited to active haploid cell promoters together with Tbp, Taf7l and RNA polymerase II. RNA-seq analysis identifies a set of genes activated in haploid spermatids during the first wave of spermatogenesis whose expression is down-regulated by Trf2 inactivation. We therefore propose that Trf2 is recruited to the preinitiation complex as a testis-specific subunit of TFIIA/ALF that cooperates with Tbp and Taf7l to promote haploid cell gene expression.